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ABSTRACT 

Long gamma-ray bursts (GRBs), among the most energetic events in the Universe, are explosions of massive and short¬ 
lived stars, so they pinpoint locations of recent star formation. However, several CRB host galaxies have recently been 
found to be deficient in molecular gas (H 2 ), believed to be the fuel of star formation. Moreover, optical spectroscopy 
of CRB afterglows implies that the molecular phase constitutes only a small fraction of the gas along the CRB line- 
of-sight. Here we report the first ever 21 cm line observations of CRB host galaxies, using the Australia Telescope 
Compact Array, implying high levels of atomic hydrogen (Hi), which suggests that the connection between atomic 
gas and star formation is stronger than previously thought, with star formation being potentially directly fuelled by 
atomic gas (or with very efficient Hl-to-H 2 conversion and rapid exhaustion of molecular gas), as has been theoretically 
shown to be possible. This can happen in low metallicity gas near the onset of star formation, because cooling of gas 
(necessary for star formation) is faster than the Hl-to-H 2 conversion. Indeed, large atomic gas reservoirs, together with 
low molecular gas masses, stellar and dust masses are consistent with CRB hosts being preferentially galaxies which 
have very recently started a star formation episode after accreting metal-poor gas from the intergalactic medium. This 
provides a natural route for forming GRBs in low-metallicity environments. The gas inflow scenario is also consistent 
with the existence of the companion H 1 object with no optical counterpart ~ 19 kpc from the CRB 060505 host, and 
with the fact that the H 1 centroids of the CRB 980425 and 060505 hosts do not coincide with optical centres of these 
galaxies, but are located close to the CRB positions. 

Key words, dust, extinction - galaxies: ISM - galaxies: star formation - gamma-ray burst: general - radio continuum: 
galaxies - radio lines: galaxies 
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1. Introduction 

Long (duration > 2 s) gamma ray-bursts (GRBs) have 
been shown to be related to the collapse of massive stars 
(e.g. iHiorth et al.ll2003t IStanek et al.ll2003li . Because of the 
very short main-sequence lifetime of such stars, GRBs are 
expected to trace galaxies with on-going star formation. 
Whether GRBs and their hosts can be used as an unbiased 


the fuel of star formation (iCarilli fc Walter] l2013lh and 
the short duration of the main-sequence phase of GRB 


tool to trace star formation in the Universe (lYiiksel et al. 

2008 

Kistler et al. 2009 

: Butler et al.l 2010t Elliott et al. 

2012 

Robertson & Ellis 

l2012f) is still a subject of de- 

bate. Some results sugge 
with the general populat 

st that GRB hosts are consistent 

ion of galaxies (|Michalowski et al. 

2012bt Hunt et al. 2014allSchadv et al.l 2014 Greiner et al. 

2015] 

iKohn et al.l 20151b 

out some biases, especially in stel- 


2011 l2015blla|{cj: 


Schulze et alJl20 lit 


In order to establish the link between GRBs and star 
formation it is necessary to understand the properties of 
the gas content in the GRB hosts, which is the fuel of star 
formation. The neutral hydrogen in the interstellar gas of 
host galaxies has been routinely detected in absorption in 
the spectra of GRB afterglows. In this case, the detected 
column density probes the portion of the neutral gas in the 
host galaxy which is distributed in one sight line in front 
of the GRB. An Hi column density larger than 10 20 cm~ 2 
(a damped Lyman a - DLA) is meas ured in more than 
70 GRB spectra. I t was clear ea r ly on dBerger et al. 2001 ; 
Fvnbo et al.ll2006b ISavaglioll2006l : iProchaska et al. 2007aHt 


2007: D’Rlia et al. 201(1 20 

4) and are at most a few per- 

cent dProchaska et al.ll2009l: 

Kriihler et al.l 20131: Friis et al.1 


progenitors implies that their hosts are currently star- 

forming (e.g. Christensen et al. 2004 Castro Cero 

n et al. 

20061 2010: Le Floc’h et al. |2006: ISavaglio et al. 

2009; 

Levesque et al. 2010a]; Svensson et al. 201(1 Kriihler et al. 

2011. 2012b| Hiorth et al.l 2012b Jakobsson et al. 

2012; 

Milvang-Jensen et al. 120121: Perlev et al.l 20131 

2015b; 


2014), and this process should be fuelled by molecular gas. 


Watson et al.l 12007 1 that these absorbers were on average 
stronger than what is typically seen in quasar spectra, also 
probing the neutral gas in high-z galaxies. It was suggested 
that this is because GRBs happen inside star-forming re¬ 
gions of the g alaxy, where the gas density is higher th an in 
a galaxy halo dFvnbo et al.1 2008; P ontzen et al. I2010T) . The 


halo contains most of the galaxies’ volume, so a sight line 
towards a distant quasar is likely to intercept it, rather than 
a dense star-forming region. The other important difference 
with high-z quasar-DLAs is an average highe r metallicity 
(|Prochaska fe Wo llel 120091 : ISavaglio et al.ll20T2lh 

Molecular hydrogen was also detected in absorp¬ 
tion in a few GRB afterglow spectra. These data sug¬ 
gest low molecular hydrogen content for GRB hosts, 
as the measured molecular gas column density frac¬ 
tions (2N (H?) / (2N (K?) + A^(Hi)) are lower than 10~ 4 


l20ll. 

However, in order to measure the total gas content 
of GRB host galaxies (as opposed to column densities 
along a single line-of-sight from GRB afterglow absorption 
spectroscopy), atomic and molecular lines need to be de- 
tected in emi s sion. A fter numerous uns u ccessful searches 
(iKohno et al.l 120051: Endo et al. l2007t lHatsukade et al.1 
120071 12011b Stanwav et al. 2011 1. the carbon monoxide 
(CO) line emission from three z ~ 0.1-0.8 GRB host 
galaxies was detected and claimed to be weak compared 
with their star format i on rates (SFRs) and dust masses 
(lHatsukade et al.l 12014 IStanwav et al.ll2015blh This may 
be surprising, because molecular gas is believed to be 


Scottish Universities Physics Alliance 


However, a large fraction of the interstellar medium 
(ISM) is in the form of neu tral h ydroge n. This is especially 
true for dwarf galaxies (e.g. lHunt et al.ll2f)14bjb and to some 
extent for moderate-mass galaxies (M* ~ 10 1 °-10 11 Mq; 
iLerov et "all 120081 : iLagos et ah 201 lib which d ominate the 
cosmic SFR density ( Brinchmann et all 1200411 , so we ex¬ 
pect them to host a large fraction of GRBs. Up to date, 
there has been no systematic study of the molecular and 
neutral gas content of GRB hosts, and of its relation to 
their (usually) moderate SFRs and dust masses. 

The objectives of this paper are to: i) present the first 
detection of the neutral hydrogen gas emission in a sample 
of GRB host galaxies; ii) test whether the amounts of gas 
are at odds with what is usually found in other star-forming 
galaxies; and in) consider what the atomic gas masses of 
GRB hosts can tell us about the fuelling of star formation 
and the evolutionary stage of GRB hosts. 

We use a cosmological model with Hn = 70 k m s -1 
Mpc -1 , Ha = 0.7, and = 0.3 and the IChabrierl (|2003f) 
initial mass function (IMF), to which all SFR and stellar 
masses were conve rted (by d ividin g by 1.8) if given origi¬ 
nally assuming the ISalpeterl (Il955l ) IMF. 


2. Samples and Data 

2.1. GRB hosts 

Our HI target sample comprises all five z < 0.12 GRB hosts 
in the southern hemisphere. This ensures that they are vis¬ 
ible by the Australia Telescope Compact Array (ATCA) 
and that the Hi line is in the accessible frequency range. 
Table |T] shows the basic properties of these hosts, whereas 
Table [2] lists the details of the ATCA observations. 

We performed radio observations with AT CA usi n g the 
Comp act Array Broad-band Backend (CA BB: 1 Wilson et al.l 
12011 1 on 12 Apr 2012, 18-25 Jul 2013 and 11-14 Apr 
2014 (project no. C2700, PI: M. Michalowski). The array 
was in the 1.5B configuration with baselines up to 1286 
m for the GRB 980425 host and in the 6A configuration 
with baselines up to 5939 m for other hosts. The data re- 
duction and analysis were done using th e Miriad package 
(ISault, fc Killeenlf2004 ISault, et aD 1 19951) . 

One intermediate frequency (IF) was centred at 2.1 GHz 
with a 2 GHz bandwidth (2048 channels 1 MHz wide each). 
We analysed the data separately in four 0.5 GHz frequency 
ranges centred at 1324, 1836, 2348, and 2860 MHz. For the 
GRB 980425 host we used a 50" diameter aperture to mea¬ 
sure the fluxes. Other hosts are not resolved, so we applied 
Gaussian fitting. The continuum data fo r the G RB 980425 
host was presented in iMichalowski et all (l2014bT) . 

A second IF was centred at the Hi line in the ATCA 
CABB ‘zoom’ mode with 32 kHz resolution. We subtracted 
the continuum to obtain the continuum-free data, and made 
the Fourier inversion to get a data cube with a velocity res- 
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Table 1. Properties of our sample of GRB hosts. 


GRB 

^opt 

Ref 

SFR IR / Rad 
(M 0 yr -1 ) 

Ref 

SFRuv 

(M 0 yr -1 ) 

Ref 

log (M„/M 0 ) 

Ref 

log(M d /M 0 ) 

Ref 

log(A7 H2 /M 0 )“ 

Ref 

Met 

Ref 

980425 

0.0085 


0.26 

_7 

0.14 

£ 

8.68 

2 

6.57 

_7 

< 8.05 

12 

8.60 


031203 

0.1050 

_2 

2.8 

8 

2.4 

10 

9.47 

10 

< 8.00 d 

_8 



8.27 

14 

060505 

0.0889 

3 

0.69 

T 

1.1 

10 

9.64 

11 





8.30 

11 

100316D 

0.0591 


1.7 

t 

0.30 6 


8.93 c 

t 





8.30 

15 

111005A 

0.0133 


0.42 

t 

0.16 6 

t 

9.68 

t 

6.57 

t 



8.70 

16 


Notes. Unless otherwise stated, the properties are adopted from the references given in the Ref columns. The metallicities are 12 + O/H, where the 
solar value is 8.69 JAsplund et al1l2004f) . ^ Converted to the CO-to-H 2 factor of aco = 5. ^ We calculated the SFR from the ultraviolet flux using 
the conversion of IKennicuttl dl998al) . ^ We calculated the stellar mass from ultraviolet, optical and near-infrared fluxes reported in IStarling et ahl 
d2Ql If) , ICano et al.l d2011|) and I Olivares E. et al.| d2012D using Grasil (see Sec. We calculated the dust mass using the submm flux limit and 

assuming Td — 30 K and 8 — 1.5. ___ ___ ___ 


References: 1: this work, 1: Tinnev et al.l (119981), 2: Prochaska et al.l (2004D, I3S: lOfek et al. (!2006l), |4|: IVergani et al.l (20101), |5|: IStarling et al. (201 ID 

ILevan et al.1 (20111), 7|:|Michalowski et al.l 

2014b 

), 

S: 

Watson et al. (20111), 1 9t 1 Michalowski et al. I 

20091), 10: Castro Ceron et al.| (120101), lit 

Thone et al.l (120081), 1121: Hatsukade et al. 

(2007 

), 

13 

: ISollerman et al.l ( 20051), [141: Levesaue et al 

(2010a]), [15l: iLevesaue et al.tl ( 20^l),ll6l: 


TMichalowsk^etar!r?2oT^ . 


Table 2. ATCA radio observation details. 


GRB 

Date 

tint" 

(hr) 



Beam size (") b 




1.38 


1.86 

2.35 

2.80 

Hi 


980425 

12 Apr 2012 

12/12 

23 x 14 

13.5 x 8.3 

12.0 x 7.5 

10.4 x 7.0 

36 x 22 

031203 

20-22 Jul 2013, 11,14 Apr 2014 

38.5/21 

8.4 x 

5.3 

6.1 x 3.9 

4.7 x 3.1 

4.0 x 2.7 

13 x 

5 

060505 

18,19,23 Jul 2013, 11-14 Apr 2014 

43.5/25 

13.4 x 

4.8 

9.4 x 3.4 

7.2 x 2.7 

6.1 x 2.4 

19 x 

6 

100316D 

22-25 Jul 2013, 12,13 Apr 2014 

40/19 

6.5 x 

5.1 

4.7 x 3.7 

3.7 x 2.9 

3.2 x 2.5 

8 x 

7 

111005A 

18,25 Jul 2013, 11-14 Apr 2014 

13.5/5.5 

29 x 

4.3 

17.7 x 3.3 

13.5 x 2.6 

11.6 x 2.2 

314 x 16 c 


Notes. U) The first number is the total integration time on source, and the second is the time when the Hi line was covered. 

At the frequen cy in GHz specified in the h eader, or at the f reque ncy of the Hi line. ^ The Hi data for the GRB 111005A host 
was adopted from iTheureau et all (119981) and lSpringob et all (120051 1. 


Table 3. The host photomet ry for the GRB 100316D host 
using the GROND data from lOlivares E. et al.l (|2012f) . 


Wavelength magAB 


(pm) 


0.45870 

18.09 

± 

0.03 

0.62198 

17.85 

± 

0.03 

0.76407 

17.75 

± 

0.03 

0.89896 

17.66 

± 

0.04 

1.23992 

17.49 

± 

0.05 

1.64684 

17.49 

± 

0.06 

2.17055 

17.78 

± 

0.09 


olution of 6.6 km s -1 . We then made a CLEAN deconvo¬ 
lution down to ~ 3cr, after which we restored the channel 
maps with a Gaussian beam. 

For the GRB111005A h ost we used arc h ival N ancay 
telescope Hi data from ITheureau et al.l (|l998l) and 
iSpringob et al.l (|2005ll . This galaxy was not recognised as 
a GRB host at the time of these observations. The galaxy 
appears point-l ike in th ese da ta, because of the beam size of 
4' x 20' (ISpringob et al.ll2005l) . The spectrum has a channel 
width of 25 kHz (~ 5 km s -1 ). 

The potential association of GRB111005A with the 
gala xy ESQ 58 0-49 at z = 0.01326 was suggested 
by iLevan et al.1 (1201 lh. and confir m ed by our multi¬ 
facili t y cam paign ( Xu et all l2011bl lat I Michalowski et all 
l201ll I2015I) . We measured its UV emission on the maps 
from the GALEX dMartin et al.ll2003l I20051P1 archive, ob¬ 
taining fluxes of 524 ± 38 /xJy and 79 ± 20 /iJy at the NUV 
and FUV filters, respectively. 


We calculated the stellar mass of the GR B 10 031 6D 
host using the p hotometry from IStarling et al.l (|2011h . 
ICano et~al~ (|201 1|1 and the photom etry from the data pre¬ 
sented in lOlivares E. et all (2012). The latter measure¬ 
ments are presente d in Table El We applied t he SEP fittin g 
method detailed in lMichalowski et~aI1 (|2008l . l200ih2010alM 
12012al l2014al. see therein a discussion of the derivation 
of galaxy properties and typical u ncertain t ies) based on 
35 000 templates in the library of llgjesias^Paramoet.jlJ 
(120071). plus additiona l temp lates of ISilva et al.1 (|l998jl 
and Michalowski et al.1 (1200811 . all developed in Grasiio 
(ISilva et al. Il998lh They are based on numerical calcula¬ 
tions of radiative transfer within a galaxy that is assumed 
to be a triaxial system with diffuse dust and dense molec¬ 
ular clouds, in which stars are born. 

For the GRB 031203 host we calculated a dust mass up¬ 
per limit from the 870 /im non-detection of IWatson et al.l 
(l201ll) . an d we did not use the v alue of 10 4 27 M 0 cal¬ 
culated by ISvmeonidis et al.l (|2014D based only on detec¬ 
tions at wavelengths shorter than 100 /an, as this reflects 
the amount of hot dust. Cold dust (dominating the total 
dust mass) does not emit significantly at these wavelengths. 
Using th e relation between c old and warm dust masses de¬ 
rived bv llzotov et al.1 (12014 their Fig. 1 2), the hot dust 
mass derived bv ISvmeonidis et all (120141) corresponds to 
the cold dust mass of 10 6 8 M 0 , an order of magnitude lower 
than our upper limit. 


Metallicities of our sample are shown in the second to 
last column of Table |T| They span the range 0.4-1.0 solar. 
They were derived from emission-line diagnostics, so they 
reflect gas-phase metallicities. 


1 Galaxy Evolution Explorer, http://galex.stsci.edu/ 


2 www.adlibitum.oat.ts.astro.it/silva 
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Fig. 1 . Integrated Hi line map of th e GRB hosts ( contours ; positive/negative values are solid/dashed at —2, —1, 2, 3 , 
4cr, ...) on top of the optical images (ISollerman et al.lf2005t iMazzali et al.ll20dfil : iThdne et al.ll2008t IStarling et alJl201lh . 
The spectrum was integrated in the frequency range shown as dotted lines in Fig. [2j North is up and east is to the left. 
The images are 30" x 30" except for GRB 980425, for which it is 100" x 100". The scale is indicated on each panel. The 
FWHM beam sizes of the HI data are shown as filled ellipses. The blue circles show the GRB positions. The map for 
the GRB 111005A host is not shown, as we use archival data with poor spatial resolution (Sec. 12.11) . To estimate the Hi 
properties for the GRB 060505 host only the central source was used. The north-western source is at the same redshift 
(Fig. [2]) and is analysed separately (Tableland Sec. 14.31) . 


Unless otherwise stated, when we refer to the SFRs of 
GRB hosts we use the infrared or radio estimates from the 
forth column of Table [Q 


2.2. Other galaxy samples 

In order to place the GRB hosts in the context of gen¬ 
eral galaxy populations we compared their properties 
with those of the following galaxy samples, chosen based 


on t he availability of the gas mass estimates: local spi¬ 
rals (jPevereux fo Younel Il990l) . optic al flux limited sp i- 
rals and irregulars with IRAS data (1 Young e t al.l I1989H 
local Luminous Infrared Galaxies (LIRGs: ISanders et al.l 


1991]^ l ocal Ultra L uminous Infrared Galaxies (ULIRGs; 


Solo mon et al.lll997T). the Herschel Reference Survey (HRS; 
BoseUiet_aL 1 120101: [Cortese et a .1 [ 201 2|. |2014t iBoselli et all 
2014 iGiesla et al.1 2014), Hi-selected z < 0 .02 galaxies 


dDovle fc Drinkwaterl 20061. Hi-rich galaxies (I Wang et al.l 
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Table 4. HI properties of GRB hosts. 


GRB 

(1) 

ZHI 

(2) 

VFWHM 

(km s -1 ) 

(3) 

F peak 

(mJy) 

(4) 

Tint 

(Jy km s -1 ) 

(5) 

log(±Hl) 

(K km s -1 pc 2 ) 
(6) 

log(Mm) 

(M®) 

(7) 

980425 

0.008607 ± 0.000005 

69 ±4 

27.6 ± 1.0 

2.17 ±0.09 

10.679 ±0.018 

8.849 ±0.018 

031203 



< 2.2 

< 0.69 

< 12.404 

< 10.613 

060505 

0.089508 ± 0.000034 

101 ± 32 

1.2 ±0.2 

0.15 ±0.03 

11.577 ±0.074 

9.780 ± 0.074 

060505comp 

0.089527 ± 0.000025 

71 ±15 

1.4 ±0.2 

0.14 ±0.02 

11.531 ±0.072 

9.734 ± 0.072 

100316D 



< 1.5 

< 0.30 

< 11.540 

< 9.731 

111005A 

0.013229 ±0.000014 

251 ± 10 

27.0 ±0.9 

7.75 ±0.25 

11.606 ±0.014 

9.778 ±0.014 


Notes. (1) GRB number. (2) Redshift determined from the Gaussian fit to the Hi spectrum. (3) Full width at half maximum 
of this Gauss ian. (4) Peak o f this Gaussian. (5) Integrated flux within 2a of the Gaussian width. (6) Hi lin e luminosity using 
equation 3 in I Solomon et all 11997 1. (7) Neutral hydrogen mass using equation 2 in lDevereux fc Younel (1990). The row marked 
‘060505comp’ corresponds to the companion object Northwest of the GRB 060505 position (Fig. [T]). 


l2013h . the GALEX Arecibo SDSS Survey (GASS) 
0 .025 < z < 0.05 galaxies with log(M*/M Q ) > 10 
(ICatinella et al.l 120101 ), the CO LD G ASS survey sup- 
plementing CO data dSaintonge et al.l 120111) The H I 
Nearby Galaxy Survey (THINGS: IWalter et al .1120081) . Hi- 
d ominated, low-m ass galaxies and large spiral galaxies 
(|Lerov et al.l 120081) . 0.01 < z < 0.0 3 mass-selected galax - 
ies with 8.5 < log(M*/Af 0 ) < 10 (|Bothwell et al.ll2014fl . 
dwarf Local Irregulars That Trace Lumin osity Extremes 
LITTLE) THIN GS dHunter et al.l l2012h . local dwarfs 
Stiln et aLll2013l). Herschel Virgo Cluster Sur vey (HeViCS; 


Davies et alJ l2010l ). Blue Compact dwarfs dGrossi et al.l 


201(1). Hi-select e d Arec ibo Legacy Fast ALFA (ALFAL FA; 


Giovanelli et al.l 120051) dwarfs dHuang ct al.l l2012h . a 
volume-limited s ample of dwarfs at d istances < 4 Mpc 
(lOtt et all 2012'; metal-poor dwarfs ( Hunt et abl l2014bL 


iLerov et al. 


20071) . meta l-poor dwarfs f rom t he Herschel 
Dwarf Galaxy Survey (ICormie r et all 2014 ), z ~ 1.5 
BzK galaxies (iDaddi et al. 20101 : Magdis et all 120111 : 
iMagnelli et al.ll2012f). a nd 1.2 < z < 4.1 submillime tre 
galaxies ( Bothwell et all [20131 : 1 Michalowski et al.ll2010afl . 

All dust mass estimates were converted to a common 
dust mass absorption coefficient Ksso/mi = 0.35 cm 2 g -1 , 
all ste llar masses and SFRs were converted to the IChabrierl 
( 120031) IMF, and all molecular masses were converted to 
otco = 5 Mq (K km s^ 1 pc 2 ) -1 . This Galactic value is ap- 
ropriate for 0,4-1 solar metallicity gal axies discussed here 
Bolatto ct al. l2013i : iHunt et aJll2014b h. 


SFR estimates of other galaxies are often derived 
from various diagnostics (UV, Ha, IR, radio), but they 
were shown to be bro adly consistent (1 Salim et al.l l2007t 
IWiiesinghe et alJ 120111) , even in dwa rf galaxies, except 
of very low S FR < O.OO1M 0 yr -1 dHuang et al.l l2012t 
iLee et all 1200 9!). not discussed here. 


3. Results 

3.1. Neutral hydrogen Hi line 

The fluxes at each frequency element were determined 
by aperture photometry with the aperture radius of 75" 
for the GRB 980425 (to encompass the entire H 1 emis¬ 
sion, see Fig. [T]) and of 8" for the remaining targets. 
For t he GRB 111005A the sp ectru m was directly availabl e 
from iTheureau et al.l (|1998T) and ISpringob et al.l (120051) . 
Gaussian functions were fitted to the spectra (Fig. [2]) and 
the parameters of the fit are reported in columns 2-4 of 


Table [I] The H 1 emission map derived from the collapsed 
cube within 2er from this fit (dotted lines on Fig. [2]) is shown 
on Fig. [T] This range was also used to obtain integrated Hi 
emission (Fi nt in Jykm s -1 ) directly from the spectra (not 
from the Gaussian fit, which in some cases is not a per¬ 
fect representation of the line shape). The line luminosity 
(L(j T in Kkm s ~ 1 pc 2 ) was calculated using equation 3 in 


ISolomon et al. (1 9971) and tr a nsform ed to Mhi using equa- 
tion 2 in lDevereux fc Youngl (|l990l) . 

We detected the H 1 emission of the hosts of 
GRB 980425, 060505 and 111005A. This is the first time 
when atomic gas is detected in emi ssion from a GRB 
host (see also lArabsalmani et al.ll2015l on the GRB 980425 
host, based partially on the same data). Moreover, for 
GRB060505 we detected Hi emission ~ 10" (19kpc at 
its redshift) north-west of the centre of the host H 1 emis¬ 
sion (R.A. = 22:07:02.9, Dec. = —27:48:46.4; source 
‘060505comp’ on Fig. [2] and in Table [I]). The redshift of 
this object is consistent with that of the GRB 060505 host 
(velocity offset ~ 6 ± 18 km s _1 ). This object is discussed 
in Sec. roi Its flux was not included in the estimate for the 
GRB 060505 host. 

The beam sizes of the radio continuum maps and Hi 
maps are shown in Table [2j Only the GRB 980425 host 
is resolved, and its dynam ical properties are presented in 
lArabsalmani et all (120151) . The lack of higher-resolution 
data does not affect our results, as we analyse the total 
Hi content, not its distribution. 

The atomic hydrogen masses of GRB hosts and other 
galaxies as a function of their SFR, stellar mass and 
dust mass are shown in Figures [31 [4] and [5l respec¬ 
tively. We find that the GRB 980425 host has a molec¬ 
ular gas mass fraction of Mh 2 /(Mh 2 + Mhi) < 14%, 
which is within the range for other star-forming galax¬ 
ies (a few to a few tens percent: lYoung et al. _ 19891: 


erov ct al.lfePOSHSaintongc et al.l 


Devereux fc Youn 


gll990l:lL 

lll20l4II 


20 111 : ICortese et al.l 1 20 14 iBoselli et al.ll2014r) . 

The Hi spectra of both GRB980425 and 111005A ex¬ 
hibit a double-peaked profile characteristic for a rotating 
disk. 

For ESQ 580 -49 (the GRB111005A host) 
ITheureau et all ( 1998 ) provided total widths of the H 1 line 
at 20% and 50% of the peak flux of 284 ± 15 km s _1 and 
272 ± 10 km s _1 , respectively, slightly higher than the 
FWHM given in Table [4j due to the fact that the Gaussian 
function does not represent the line profile accurately. Our 
estimate of the integrated flux (which does not involve 
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assumptions on the line shape) , agre es with 7.6 ± 0.9 Jy 


km s 1 given bv iTheureau et al~ (1998 ). 


3.2. Radio continuum 

The radio images at all four frequencies are shown in 
Fig. [6] All hosts are detected, but the significance for 
the GRB 060505 host is only around 3cr. For each host 
we fitted a power law to these fluxes obtaining the radio 
slope a and the rest-frame 1.4 GHz flux, which__was trans¬ 
formed to S FRraHio using the conversi on of iBelll (120031) . 
Finally, as in iMichalowski et, al.l (l2012blh we calculated an 
approximate measure of the ultraviolet (UV) dust atten¬ 
uation Auv = 2.51og(SFR ra dio/SFRuv), which we con¬ 
verted to visual attenuation assumin g an SMC extin ction 
curve, which gives Ay = Ayy/2.2 (|Gordon et al.ll2003h . All 
these estimates are reported in Table [5] Due to a similar 
shape of all standard extinction and attenuation curves be¬ 
tween 0.275 fi m and U-band, our Ay estimate would only 
be ~ 10% higher if we used the LMC or the Milky Way 
extinctio n curves iGordon etal.1 [20031) . and ~ 20% higher 
with the lCalzetti et ahP 2000( 1 "attenuation curve. 

Up to date 89 long-GRB hos ts have bee n tar- 


getec 

2003 

1 in tne radio con 

IVreeswiik et al.l 12001 

;miium merger et aui zuujlil 

Fox et al. 2003t Frail et al. 

2003 

: van der Horst et al. 

20051 IWiersema et al.l 

2008; 

Michalowski et al. 2009, 20124 Stanwav et al. 201(A 

2014; 

Watson et al. 20111: Perlev & Perlev 20131: Perlev et al. 

2015bl). and only 15 (~ 

17%) were detected 

the 

hosts of GRB 980425 (IMichalowski et al. 20091). 980703 

(iBerger et al.) 2001|), 000418 

. 010222. (Berger et all 20031). 

021211 (Michalowski et al. 2012bl). 031203 (Stanwav et al. 

201X1 IWatson et al. 2 Oil], 

IMichalowski et al.l 2012b!). 


(Perlev et al. 2015 


051022, 080207, and 090404 (jF 

051006, 060814, 06 112 1, 070 306, ,___ 

0805 17 dStanwav et al.l l2015al) . 100621A ( Stanwav et al.l 
l2014fl : with the addi tion of the hosts of short (dur ation 
< 2 s) G RBs 071227 (iNicuesa Guelbenzu et al]l2014l) , and 
120804A dBerger et all 2013FT 

We provide three more detections (060505, which has 
been targeted before but not detected, 100316D and 
111005A) bringing the fraction of detected hosts to 18/91 
(~ 20 %). 


4. Discussion 

4.1. Large H I reservoirs in GRB hosts: early stages of star 
formation 

It is still an unanswered question whether GRB hosts 
are consistent with the general population of star-forming 
galaxies. As shown in Figs. [3l []] and El GRB hosts at 
z < 0.12 contain large atomic gas reservoirs, within the 
ranges expected for galaxies with similar SFRs, stellar 
masses and dust masses. 

We performed a two-dimensional Kolmogorov-Smirnov 
test using th e pub licly available IDL procedures^ of 
lYoachim et alJ (j2009l) . We com pared t he GRB h ost sam - 
ple to the Hl-selected galaxies (iDovlc fc Drinkwatci 


on the SFR-Mhi plane, and the HRS (iBoselli et ah 


2006) 


2010) 


galaxies on the M*-Mhi plane, because GRB hosts span 


similar SFR and M* ranges to those of these galaxies, re¬ 
spectively. For GRB 031203 and 100316D with no Hi detec¬ 
tions, we assumed 20 different values of their Mhi between 
M*/100 (a conservative lower limit) and the derived Mhi 
upper limits (TablcEJ). We obtained p- values of 11-79% for 
the SFR-Mhi plane, and 0.2-24% for the M*-Mhi plane. 
This implies that GRB hosts have consistent (or somewhat 
higher) atomic gas masses when being compared with other 
samples. If the Hi masses of non-detected hosts are close 
to the derived upper limits, then the probability for the 
M*-Mhi plane is 0.2%, corresponding to a ~ 4tr discrep¬ 
ancy between the GRB host and the HRS samples. In such 
a case, GRB hosts are all located in the high-Mni end of 
the distribution or have low M*, which is also evident on 
Fig. [I] If these galaxies are well below the derived upper 
limits, then the GRB host population is consistent with the 
comparison sample. 

The mass at which the contribution of local galaxies to 
the Hi mass function is the highest is l og (Mhi /MM ~ 9.6 , 
as derived from the Mhi function of IZwaan et al.l (2 0051 ). 
This is similar to what we found for GRB hosts. This is con¬ 
sistent with numerical simulations, suggesting that galaxies 
with low SFRs (similar to low -0 GRB hosts) dominate the 
Hi mass densit y at all redshifts, including the local uni¬ 
verse (Fig. 6 of lLagos et al.l [20141) . In particular, galaxies 
with SFR=0.1-5Mg yr _1 and 0.1-1 M 0 yr _1 contribute 
~ 45% and ~ 30% to the Hi mass density at z < 1. 
(lLagos et al.l 12012L 12014 private communication) All our 
Hi-detected GRB hosts are in this SFR range, so they can 
be considered typical galaxies at low -0 with regards to their 
Hi content. 

These considerations indicate that, at least at 0 < 
0.1, GRB hosts can be regarded as normal star-forming 
gal axies. This is cons is tent wi th the conclus i on pres ented 
in IMich alowski et, al.l (l2012bl) . iHunt, et ~aI1 (|2014al) and 
iKohn et alJ ( 20151) based on radio and far-infrared contin¬ 
uum data. 

To compl ete th is picture we note that iPerlev et al.l 
(120131 Mm and IVergani et alJ (|2014h found that GRB 
hosts at 0 < 1.5 have lower stellar masses than what 
would be expected from the assumption that GRBs trace 
the cosmic s tar formation acti v ity in an un biased way. 
Additionally, IMichalowski et all (l2014bl . 120151) found dust 
masses of two low -0 GRB hosts to be close to the lower 
envelope of other galaxies. 

Normal SFRs and atomic gas masses, together with low 
stellar and molecular gas masses are consistent with GRB 
hosts being preferentially galaxies which have very recently 
started a star formation episode. In this scenario, GRB 
hosts have not had enough time yet to use their atomic 
gas reservoir (which is then high) and to produce stars and 
dust (which catalyses the molecular gas formation). 

All our targets are in the nearby Universe, so are likely 
different from the high -0 counterparts. Only with the ad¬ 
vent of the Square Kilometre Array (SKA) the Hi emission 
of GRB hosts (and other galaxies) at higher redshifts can 
be detected. 


4.2. H 2 vs Hi: what is the fuel of star formation? 

Large atomic gas content, as measured from our data, to¬ 
gether with the molecular gas deficiency can in principle be 
explained in four different ways, which we discuss below, 


http://www.astro.Washington.edu/users/yoachim/code .phjproviding evidence that the first two contradict observa- 
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Table 5. Radio continuum properties of GRB hosts. 


GRB 

±1.38 

(/Uy) 

±1.86 

(pJy) 

± 2.35 

(tUy) 

± 2.80 

(pJy) 

-^restl.4 

(/^Jy) 

a 5 

SFR mdio c 
(M @ yr^ 1 ) 

Ay d 

(mag) 

980425 

840 ± 160 

1110 ±210 

710 ± 190 

400 ± 200 

912 ± 263 

-0.58 ±0.41 

0.11 ±0.03 

~ 0 

031203 

383 ± 20 

324 ± 11 

319 ± 12 

324 ± 13 

363 ± 26 

-0.19 ±0.09 

2.83 ±0.20 

0.08 

060505 

76 ± 35 

67 ± 18 

45 ± 18 

47 ± 15 

85 ± 50 

-0.81 ±0.72 

0.69 ± 0.40 

~ 0 

100316D 

657 ± 21 

531 ± 13 

441 ± 13 

403 ± 13 

681 ± 33 

-0.71 ±0.06 

1.73 ±0.08 

0.86 

111005A 

245 ± 30 

192 ± 18 

160 ± 16 

124 ± 16 

252 ± 46 

-0.92 ±0.22 

0.08 ±0.02 

~ 0 


Notes. Flux densities are given at the observed frequencies indicated in GHz in the column header. U) Rest-frame 1.4 GHz 
flux density from the power-law fit to the data. ^ Radio power-law slope. G) Using the conversion of iBelll (2003). G) Visual dust 
attenuation calculate d from the ultraviole t attenuation Auv = 2.5 log(SFR ra di 0 /SFRuv) assuming an SMC extinction curve, which 
gives Ay = Any/2.2 (iGordon et al.ll200 3). Ay ~ 0 mag is reported if a formal value of SFRuv exceeds the SFR. ra di 0 . 


tions and that the last one (star formation directly fuelled 
by a recent accretion of atomic gas) is the most plausible 
(because of the metallicity considerations and the success 
to explain low stellar masses of GRB hosts). 

First, it is possible that very strong UV radiation 
from the star-forming region in which a GRB progeni- 
tor is born dissipated the available molecula r hydrogen 
(lHatsukade et al.l 12014 IStanwav et al.1 l2015bl) . However, 
this is unlikely to affect the entire galaxy, but just the 
vicinity of the GRB. Moreover, such process would also de¬ 
stroy dust, whereas the molecule-poor regions c lose to the 
GRB sites are dust-rich (lHatsuka de et al.l 2014fl. as is th e 
region close to the GRB 980425 ((Michalowski et ahll2014bf) . 
Finally, very intense UV radiation dissipating the remain¬ 
ing molecular gas would imply that the surrounding of a 
GRB is at the end of the intense star-formation episode, 
which should result in an enha nced metallicity in that 
region, contrary to observations dGhristensen et al.l 120081: 


Modiaz et al . 2008tlThdnc et al.ll2008ll20l4lLevesaue et al.l 


201 Oal 120111 lHan et al.l 12010 1. 

Second, low molecular gas masses in GRB hosts were 
derived from the CO lines, which may not be a good 
tracer of molecular gas at low metallicities if dust shielding 
is weak dBolatto et al.l [2013h . and if a significant amount 
of H 2 gas is CO-dark, as is the case for metal-poor 
dwarfs. However, this is not the case for GRB hosts for 
which low molecular gas content was claimed, because they 
have 0.5-1.0 solar metallicity (ICastro-Tirado et al.1 20071 
Graham et al.j | 2004 iLevesoue et al.l 2010bl Stanwav et al.l 
2015al1 . Similarly, in our sample the metallicity is not ver 


low either, around 0 . 4-1.0 solar ( Sollcrman^_alJ_|200 


Levesque et al.l l2010al iThonc ct~al 


20081 Levesque et al l 


20nT see Tabled]). Moreover, a low molecular gas content 


dFvnbo et al. 2006 

: Tumlinson et al.l 20071; Kriihler et al. 

201.4 D’FJia, et al l 

20 lj- This is not due to molecule de- 


below), so this explanation is possible only if some prop¬ 
erties of GRB hosts make the H l-to-H 2 conversion unusu- 
ally effective, w hich is not accounted for in the models of 
iKrumhoQ ( 2012 ). 

Finally, the fourth and our preferred option is that star 
formation is directly fuelled by atomic gas from the inter- 
gala ctic medium. This was theoretically sh own to be possi¬ 
ble dGlover &; Clarfdl2012t lKrumhoT3l2012f) , and must have 
been the case for the very first stars in the Universe (in the 
absence of dust and molecular gas). This can also take place 
at any redshift, even in a galaxy not particularly metal-poor 
on average, as long as it happens at early stages of a star 
formation episode, when the first stars are born in a col¬ 
lapsing cloud formed out of newly accreted metal-poor gas. 

Physically, this scenario rests on two points. First, 
observations indicate that star formation begins within 
roughly one free-fall time of the appearance of a cold, 
d ense phase of th e ISM (for example see the review by 
iDobbs et al.ll2014l) . Only the duration of star formation is 
a matter of a debate, but this is not relevant here, as we 
concentrate on the beginning of a star formation episode. 
Indeed, some observations and theoretical models indicate 
that star formatio n finishes quickly after it starts, around 
1-2 free-fall times (lElmegreenl 200C; Hartmaim^jd. 2001 


Jeffries et al.l 1201 it Reggiani et al.1 20111 : Dobbs fe Pringlel 

201311 . whereas other suggest that the duration of star for- 


I 


mation can be significantly longer (several up to 10 free- 
fall times), though it still starts around one free-fall time 
jTme^iLL 2006 • Krumholz et al.l 12006 ; iKawamura et al.l 
120091: iFukui ei~aLl(l999l1.~ 


in GRB hosts has also been found independently of CO ob¬ 
servations, using the optical afterglow spectroscopy, tracing 
molecular gas column density along the GRB line-of-sight 


struction by the GRB itself, as thi s can onl y be effect ive up 
to a few parsecs from the GRB (jDraine fc Haoll2002f) . 

Third, atomic gas may be converted into the molecular 
phase and be immediately used for star formation (which 
explains the lack of molecular gas). However, the timescale 
of the H l-to-H 2 conversion is usually longer tha n the cool¬ 
ing t imescale (and the star formation timescale; lKrumhoI3 
12012 ). so this could happen only at approximately solar 
metallically, at which Hi-to-H 2 conversion is quicker than 
the star formation timescale. Metallicities of GRB hosts 
and GRB sites are usually much lower (see the discussion 


Second, the timescale for the appearance of a cold phase 
is determined by the gas cooling time, while the appearance 
of H 2 is determined by the chemical equilibration time. At 
low metallicities the cooling time is longer, but still re¬ 
mains much shorter than the free-fall time, w hereas the 
chem ical equilibration time can be much longer (iKrumholzl 
120121 . fig. 1), leading to the formation of a cold phase and 
the onset of star formation before gas is able to convert 
from Hi to H 2 . Specifically, theoretical models (IKrumholzl 
12012 ) predict that at the metallicities of ~ 0.3 solar typ- 


Christensen et al.ll2008l 

Modiaz et al.l 2008k Levesaue et al.l 

2010al. 201 it Han et al. 

201(4 Thone et all 2014ll gas can 


cool to low temperatures in less than a free-fall time, but 
that full conversion to molecular gas requires close to ten 
free-fall times, so the first stages of star formation are 
fuelled by atomic gas. Moreover, the observed GRB site 
metallicities could even overestimate the metallicity when 
the GRB progenitor formed, since the metallicity will have 
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been increased during the progenitors lifetime by on-going 
mixing with the more metal-rich gas existing in the hosts. 
If so, the discrepancy between molecule-formation and free- 
fall timescales would be even larger. Thus our preferred sce¬ 
nario is that GRB progenitors potentially form in the first 
burst of star formation that takes place in newly accreted 
low-metallicity inter galactic medium (IGM) gas. The pro¬ 
genitor then would end its life as a GRB before there is 
time for substantial H 2 to form, explaining the paucity of 
H 2 absorption features seen in GRB afterglows, and the 
paucity of CO emission from the host galaxy. 

Indeed recent star-formation and large Hi content for 
other galaxies (as also observed for GRB hosts) are be¬ 
lieved to be a result of recen t inflow of met al-poor gas 
from the intergalactic medium (iDave et al.ll2013h . This cre¬ 
ates metal- poor regions in a galaxy (as observed in other 
galax ies bv lGresci et al.ll2010l : ISanchez Almeida et al.ll2013l 
l2014f) , which mix with the surrounding metal-rich gas over 
a relatively long timesca l e (< 100 Mvr: lYang fc Krumholj 
l2012t iGrand et al.l 120151 : iPetit et alll2015f) . Indeed. GRBs 
were shown t o explode in the most metal-poor regions 
of their hosts dCh ristense n et al.1 [20081 : IThone et al.1120081 
l20l4lLevesaue et al.l201 it) . The GRB site metallicities are 
~ 0.3 solar, which would imply that less than half of the 
accreted gas converts to t he molecular phase after one free- 
fall time (lKrumhoI32012f) . This fraction can be lower if the 
accreted gas had in fact lower metallicity than the GRB 
site measurements, which have been increased by on-going 
mixing with the more metal-rich gas existing in the hosts. 
We note that the bulk of the star formation in these regions 
(at later stages) will have to be fuelled by molecular gas, 
when the conversion from the atomic to molecular phase is 
complete. 

The recent accretion of ato mic gas also explains hig: h 
specific SFRs of GRB hosts ie.g. lGastro Geron et al.l[201Qlh 
because the enhanced SFR has not lasted very long, so the 
stellar mass of a galaxy is lower than what would be ex¬ 
pected from its high SFR. Moreover, this mechanism di¬ 
rectly predicts that the regions around GRBs (at which 
the gas is accreted) have enhanced SFRs and dust masses 
compared with other parts of the hosts, consistent with ob¬ 
serva ti ons dLe Floc’h et al.l |2006l |2012t IChristensen et al.| 


20081 iThone et ah 2 00S . 2014 Levesque et alJ 20lit 


Hatsukade et al. 2014 : Michalowski et al.l 2014blh 

Our data therefore is consistent with a scenario whereby 
GRBs are preferentially produced when low-metallicity gas 
accretes onto a galaxy and undergoes rapid cooling and 
star formation before it either forms H 2 or mixes with the 
higher-metallicity gas in the remainder of the galaxy. This 
mechanism can be tested with positional and velocity in¬ 
formation from high-resolution Hi observations, which may 
reveal concentration of HI close to metal-poor regions. 

If confirmed, this mechanism will provide a natural ex¬ 
planation of the low-metallicity and 1ow-Mh 2 preference in 
the framework of the GRB collapsar model, which requires 
that GRB progenitors have low metallicity in order to re¬ 
duce the loss of m ass and angular momentum (required for 
launching the .jet; Yoon fc Langed 120051 : lYoon et ahl 120061 : 
IWooslev fc Hegerll2006lh In contrast, at later stages of star 
formation molecular gas is the dominant phase in the in¬ 
terstellar medium, but the metals are well mixed, and gas 
has been further enriched, so massive stars do not end their 
lives as GRBs, and such metal- and molecular-rich galaxies 
do not become GRB hosts. 


However, this picture is complicated by the fact 
that some GRBs have been found in ga laxies with 


solar or supersolar metallicities dProchaska et a~ 20091 


Levesaue et al.1 20 

LObt Kriihler et al.l 2012a 

Savaglio et al. 

2012): 

Elliott et al.l 

2013! ISchulzc et al.l 2014 

Stanwav et al. 

2015a 

; Hashimoto et al|2015|; Schadv et al. 

20151). This can 


be explained in two ways. First, these metallicity measure¬ 
ments were obtained for the entire galaxies, so do not rule 
out local metallicity decrements (as observed for other GRB 
hosts), predicted by the mechanism of recent metal-poor 
gas accretion. 

Second, this high-metallicity problem is alleviated by 
some models, which predict the GRB preference for low- 
metallicity, but do not exclude metal-rich examples, for ex¬ 
ample the model in which a GRB progenitor is a rapidly 
rotating Helium star created by a merger of post-mai n 
sequence stars (Ivan den Heuvel fc Portegies Zwartl I2Q13T ). 
In both collapsar and binary scenarios the low-metallicity 
preference (even if it is not strict) is consistent with our 
interpretation that GRB hosts have recently started a star 
formation episode. 

We note that the IMF in met al-poor envi r onments has 
been suggested to be t op-heav v dBatel 2005 Zhang et al.l 


120071 : iMarks et al.ll2012t but see iMvers et al.1 201 ID . Hence. 


the sites of GRB explosions are promising places to look for 
the top-heavy IMF. This IMF would imply higher massive 
star (and hence GRB) production per unit star formation 
rate in such environments. 


Summarising, our results provide the first observational 
indication of star formation fuelled by atomic gas. The link 
between the star formation and the atomic gas for all galax¬ 
ies is visible in Fig. [8] where the Hi mass is strongly corre¬ 
lated with SFR. The linear fit to all datapoints gives 


log(M ffl /M 0 ) = (0.76 ± 0.01) x log(SFR/M 0 yr” 1 ) 

+ (9.64 ±0.01) (1) 


The Spearmans rank correlation coefficient for this relation 
is 0.77, indicating a significant (~ 30cr) correlation. The 
scatter around this correlation is ~ 0.38 dex, which is likely 
a consequence of galaxies with similar SFRs being at differ¬ 
ent evolutionary stages. The correlation with similar slope 
was also noted using the surface densities of SFR and Mht 
dKennicuttl Il998bl : IBigiel et~ahl 120101 : iRovchowdhurv et al.l 
120141b 

The slope of ~ 0.76 (or ~ 1.3 with inverted axes) 
is shallower than the linear (slope of unity) Schmidt- 
Kenni cutt relation betw een molecular gas and SFR den¬ 
sities (iBigiel et al.ll2008 '). This indicates that overall molec¬ 
ular gas is better correlated with star formation, so is 
its major fuel. However, in outer parts of spiral galax¬ 
ies, and in Hi-dominated galaxies, the relat ion between 
SFR and H 1 column density is in fact linear (IBigiel et al.l 
l2008i- j20iQ) and the SFR-gas correlation improves when 
ato mic gas is taken i nto accou nt, not just the molecular 
gas (iFumagalli fc Gavazzill2008l their fig. 2). This was in¬ 
terpreted as Hi fuelling of low-levels of star formation in 
the outer parts of disks. A similar mechanism may also 
operate in some star-forming regions of GRB hosts, be¬ 
cause these galaxies have low molecular gas content, and 
are likely observed at the beginning of a star formation 
episode, just after accreting atomic gas from the IGM, un¬ 
like other galaxies. 
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4.3. Location of the Hi emission 

Despite the low resolution of our HI data we can infer 
useful information about the location of the Hi emission, 
which supports our gas inflow scenario. First, the compan¬ 
ion Hi object (not a nearby galaxy, as we show below) is 
only ~ 19kpc from the GRB 060505 host (Fig. Q]and[2]), so 
they must be physically relat ed ( the optical rad ius of the 
GRB 060505 host is ~ 11 kpc: lThdne et al.l[2008 b With the 
current resolution we cannot definitely determine the na¬ 
ture of this relation, but its existence is consistent with HI 
gas inflowing from the north-western direction, either in a 
form of a star-free atomic gas cloud, or with extremely low 
stellar surface density (we discuss below that the astrom¬ 
etry of the radio and optical images are consistent). The 
radial velocities of the GRB 060505 host and this cloud are 
identical (offset of ~ 6± 18km s _1 , TablcQ]), which is con¬ 
sistent with the scenario that the cloud is flowing directly 
toward the host, or that the orbit is perpendicular to the 
line-of-sight. The existence of this objects may indicate that 
the host of GRB 060505 is at the earliest stage of the in¬ 
flow compared with other hosts, so the inflowing material 
is still visible outside the host, unlike for other members of 
our sample. The inflow scenario for this cloud needs to be 
tested by high-resolution HI observation, which can reveal 
its morphology and velocity structure. 

There are two other explanations of this feature, but 
they are unlikely. First it could be a signature of a major 
merge® with another galaxy, but the HI companion ob¬ 
ject does not have any optical counterpart (despite high 
log(MHi/M Q ) ~ 9.7; Tablc[4]), so cannot be understood as 
a companion galaxy, but rather a large gas cloud. Namely, 
it is not detected with the HST B- band image dow n to a 
(point-source) limit of 27.1 mag AB (|Ofek et al.ll2007ll . This 
corresponds to an absolute magnitude of Mb > —10.93 
mag, and a luminosity Lg < 2.9 x 10 6 Lq. Hence it has 
log(MHi/Ts) > 3.3, much higher than those of dwarf galax- 
ies (from —1.5 to 0.75; Fig. 1 of lHunter et al.ir2012H . This 
limit assumes that this object is point-like in the HST reso¬ 
lution, but it would need to be larger than the HST resolu¬ 
tion ~ 20 times (~ 2", or ~ 3kpc) to make its Mm/L b 
ratio consistent with those of dwarf galaxies. But then 
it would be an extremely low-surface brightness object. 
Moreover, no e mission l ines were fo und at this position in 
the IFU data of IThone et al.l (120141) . 

Second, the HI companion may be ejected from the 
GRB 060505 host via tidal interaction or ram pressure strip¬ 
ping. However, this explanation is unlikely given the lack of 
nearby companions and that it is not located in the high- 
density cluster environment ISec. fPl) . 

Another information supporting our gas inflow interpre¬ 
tation is that for both the GRB 980425 and 060505 hosts 
the centre of the HI emission is shifted away from the opti¬ 
cal centre of the galaxy towards the GRB position (Fig.[T|). 
If a galaxy has a steady Hi disk, its center should coincide 
with the optical center of this galaxy. On the other hand, if 
a galaxy has received a significant gas inflow from one di¬ 
rection, then in low-resolution data the HI centroid will be 
offset from the optical center towards this direction. Hence, 
the offsets for GRB 980425 and 060505 host are consistent 
with an inflow of atomic gas on the regions around the GRB 


4 This could not be a minor merger because these objects have 
similar Hi masses. 


positions. This will be verified with future higher-resolution 
observations. 

To test the ATCA astrometry we imaged the phase cali¬ 
brator for the GRB 060505 observations (2149-287) obtain¬ 
ing the position 21:52:03.735, —28:28:28.256, which is only 
~ 0.04" away from the catalogue position 21:52:03.7352, 
-28:28:28.21® We also searched for common objects in 
the radio and optical images. There is no other significant 
object in the integrated Hi map (Fig. [[]), which would 
be inside the optical map, but in the radio continuum 
image ~ V south-west of the GRB position there is an 
~ 0.4 mJy source at the position 22:07:00.093, —27:49:14.89 
with an optical counterpart at the position 22:07:00.133, 
—27:49:15.32. Hence, the radio map is shifted only 0.7" to 
the West with respect to the optical map. 

The astrometry of the optical image (IThone et al .1120081) 
was tied to the Two Micron All Sky S urvey catalogue 
dJarrett et alJl2000t ISkrutskie et al.ll2006D . and is accurate 
to a fraction of arcsec (~ 0.1"). Moreover, it is consi stent 
with more recent astrometry from lHiorth et all (2012). 

Hence we conclude that the astrometry of the optical 
and radio images are consistent, and therefore the HI- 
optical offsets for the GRB 980425 and 060505 hosts are 
real. 

Finally, we note that the con centration of atomic gas 
in the south-eastern (blue-shifted; iMichalowski et akll2015l) 
part of the GRB 111005A host is also hinted by its asym¬ 
metric Hi line profile (Fig. [2]). The blue-shifted line part is 
more pronounced, so the south-eastern part contains more 
atomic gas. The line profile for the GRB 980425 is more 
symmetric, but this is expected, because its metal poor re¬ 
gion (close to the GRB site) has the velocity close to the 
systemic one dGhristensen et al]l2008 ). 


If GRB hosts were located in dense large-scale regions 
(e. g. close to the cores of massive groups or clusters), 
then our interpretation of recent metal-poor gas inflow 
would be difficult to advocate, because of scarcity of metal- 
poor gas in such environments. Howev er, the GRB 98042 5 
host has been shown to be isolated dFolev et al.l 120061) . 
While the GRB 060505 host lies in the foreground of a fil- 
ame ntary structure a few Mpc away from a galaxy clus¬ 
ter (IThone et al.l 120081 ), at such distances gas stripping or 
evaporatio n does not take p lace, and the SFR is actually 
enhanced dPorter et al.ll2008h . 

In order to investigate the environments of other hosts 
in our sample we searched the NASA/IPAC Extragalactic 
Database (NED) for companion galaxies within 1 Mpc 
in projection and within ±1500 km s _1 in redshift. The 
GRB 031203 and 100316D hosts have no nearby galaxies 
within these search criteria. The GRB 111005A host is not 
in any NED group catalog either, and there are only 6 
nearby galaxies nearby, all > 500 kpc away. 

Hence, the host galaxies are all relatively isolated in 
terms of large-scale environment, which makes it possible 
that they experienced inflows of metal-poor intergalactic 
gas. 


5 http://www.narrabri.atnf.csiro.au/calibrators/calibrator_d 
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4.5. Little dust-obscured star formation 

Only for the GRB 100316D host the radio SFR exceeds sig¬ 
nificantly the UV-derived value, implying Ay ~ 0.86 mag 
(Tables |T] and 0. For the remaining hosts the comparable 
SFR radio and SFRyv indicate that there is very little dust- 
obscured star formation in these hosts. This is consistent 
with our interpretation that GRB hosts are at the begin¬ 
ning of the star-formation episode, when the dust content 
has not had time to accumulate to make star-forming re¬ 
gions optically thick. 


4.6. Radio spectral slope 

The radio spectral slopes (Tabic 0 of all but one hosts are 
consistent with those of star-f orming galaxies both local 
and a t high reds hifts (~ — 0.75: [Condonlll99^lDunne et al.1 

l2009t llbar etTDIiZOTnjb _ 

As suggested in iMichalowski et all (2012b), the spec¬ 
tral slope of the GRB 031203 host is significantly flat¬ 
ter (~ —0.19 ± 0.09), inconsistent with the normal value 
for star-forming galaxies. This implies a sig nificant free- 
free emission (or synchrotron self-absorption; Condonlll992l) 
and, hence, a younger st e llar population dBressan et al l 
2002t Gannon fc Skillmanl 12004 iHirashita fc Hunt! 20061 : 
Clemens et al.l l2008lb because only massive stars emit 


Lyman continuum photons, which create H II regions, re¬ 
sponsible for free-free emission. The flat slope for the 
GRB 031203 host also rules out any significant AGN contri¬ 
bution to the radio flux. This is consistent with its location 
in the star-forming part of the Bald win-Phillips-Terlevic h 
(BPT; [Baldwin et al.lll98lh diagram (I Watson et al.ll20TTlh 


5. Conclusions 

We report the results of the first ever Hi survey of GRB 
host galaxies, detecting three out of five targets. Large in¬ 
ferred atomic gas masses, together with low molecular gas, 
stellar and dust masses are consistent with GRB hosts be¬ 
ing at the beginning of a star formation episode, after ac¬ 
creting metal-poor gas from the intergalactic medium. This 
star formation may potentially be directly fuelled by atomic 
gas (or with very efficient H i-to-H 2 conversion and rapid ex¬ 
haustion of molecular gas), which can happen in low metal- 
licity gas near the onset of star formation, because cooling 
of gas (necessary for star formation) is faster than the HI- 
to-H 2 conversion. This provides a natural route for form¬ 
ing GRBs in low-metallicity environments. The gas inflow 
scenario is consistent with the existence of a companion 
Hi object with no optical counterpart ~ 19kpc from the 
GRB 060505 host, and with the offset towards the GRB po¬ 
sitions of the Hi centroids for the GRB 980425 and 060505 
hosts away from the optical centres of these galaxies. 
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Fig. 2. Hi spectra of the GRB hosts ( histograms ) with the Gaussian fit overplotted (red lines). These fits are presented 
only for illustration, as the total fluxes were determined from the direct integration of the spectra (Sec. 13.111 . Vertical 
dotted lines show the frequency range used to measure the total HI emission and to produce the integrated HI maps 
(Fig. ED- The spectrum marked ‘060505comp’ corresponds to the companion object north-west of the GRB 060505 position 
(Fig. ED. 
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Fig. 3. Hi mass as a function of infrared/radio-based star formation rate (SFR) of GRB hosts (red circles ) and other 
galaxies, as indicated in the legend and described in Sect. 12.21 The solid black line is a linear fit to all the data (eq. |T|). 
GRB hosts are consistent with the general star-forming galaxy population. 
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Fig.4. Hi mass as a function of stellar mass of GRB hosts (red circles) and other galaxies, as indicated in the legend 
and described in Sect. 12.21 GRB hosts are consistent with the general star-forming galaxy population. 
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Fig. 5. Hi mass as a function of dust mass of GRB hosts (red circles) and other galaxies, as indicated in the legend and 
described in Sect. 12.21 Dotted lines denote constant gas-to-dust mass ratios with values indicated on the left. GRB hosts 
are consistent with the general star-forming galaxy population, especially with the Herschel Reference Survey (HRS) 
galaxies with similarly measured dust masses. 
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Fig. 6. Optical (first column; ISollerman et al.ll2005t iMazzali et al.ll20d5 lTh5ne et al.ll2008t IStarling et al.ll20TTl and this 
work) and radio continuum images centred at 2800, 2340, 1860 and 1340 GHz. North is up and east is to the left. Each 
panel is 30" x 30" The FWHM beam size is shown on each radio panel as an ellipse. The blue circles show the GRB 
positions. 
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